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We study the phase structure of the massive one flavour lattice Schwinger model on the basis of the finite size 
scaling behaviour of the partition function zeroes. At /? = we observe and discuss a possible discrepancy with 
results obtained by a different method. 



1. INTRODUCTION 

Since its introduction by Schwinger in 1962 [1], 
two dimensional QED has been continually ex- 
amined and has attracted both high energy and 
mathematical physicists. The model shows im- 
portant aspects of realistic models, is exactly sol- 
vable and the Wightman axioms hold for the bo- 
sonized fields [2]. Among the properties are brea- 
king of chiral symmetry (mass generation) based 
on confinement properties, charge shielding and 
the topologically interesting fact of a continuum 
of degenerated vacua labeled by an angle 6. Cer- 
tain properties extend to the massive version of 
the Schwinger model, where chiral symmetry is 
explicitly broken [3] . The lattice version of the 
Schwinger model has served extensively as a la- 
boratory to test fermionic algorithms, mainly for 
Kogut-Susskind fermions. Its formal continuum 
limit from a constructive field theory point of 
view was examined in [4] . 

In particular the Schwinger model with fermi- 
ons on the Wilson representation may teach us 
how to approach continuum physics due to its 
analogy to lattice QCD with Wilson fermions. 
Both models are believed to have similar zero tem- 
perature phase structure: For any gauge coupling 
(3 there exists a critical hopping parameter ndP)- 
Thus continuum physics could be approached 
along this line which ends at Kc(oo) = l/(2d), 
a critical point with a second order transition. 

Up to now it is (from a rigorous point of view) 
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unproven and thus a matter of believe how far 
this line really extends into the plane. It is 
assumed and supported by strong coupling data, 
that there is a critical point at /3 = cxd with < 
Kc < oo. In the strong coupling limit expectation 
values of operators are expanded around (k, (3) = 
0. 

j dij{U)0{U)detM{K,U)e-^°'-fi'^) 

= J2n,m Cn.mK"/?™ • (1) 

Here M denotes the lattice Dirac operator. This 
expansion is known to converge at least in a non- 
vanishing domain around the expansion point, i.e. 
|k| < Pk(P) and (3 < pp [5] and has in that range 
a finite mass gap. Usually Pk(P) is then identified 
with a critical value Kc(/3) G M"*". This, however, 
is not necessarily true, since the singularity could 
occur for a complex k on the circle \k\ = pPP). 

For the Schwinger model it can be shown that 
/Ok(0) < I [6]. Further checks, e.g. numerical 
simulations, have to be done to verify the as- 
sumption Kc = Pk.- Typically one would calculate 
the chiral condensate {ipip), but for Wilson fermi- 
ons this is not an order parameter and no simple 
parameter is known. Another standard possibi- 
lity is to determine the correlation length (inverse 
mass gap) and look for divergences in the plane 
of coupling constants. Doing this in an analy- 
tic calculation (hopping expansion) the necessary 
approximations are too rough for a trustworthy 
verification whether there is a critical point. 

A direct numerical simulation of the one fia- 
vour Schwinger model for general (3 is not possi- 
ble with present means since detM(K, U) as well 
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as quadratic forms ((j>, M~^(K,U)(j>) are not po- 
sitive definite for certain gauge field configurati- 
ons. In this case Metropolis type algorithms will 
fail. It has been shown in [7] that the one flavour 
Schwinger model at /3 = can be mapped onto 
a "seven- vertex" model (the eight-vertex model 
for a particular choice of weight parameters). In 
the limit k ^ co that model becomes the criti- 
cal six- vertex model [8] . Therefore the Schwinger 
model has at least one critical point at /3 = and 
Kc(0) = cxD, far beyond the range of the cluster 
expansion. This map further offers the oppor- 
tunity for a direct simulation of the one flavour 
Schwinger model at least at /3 = and allows 
for a test of strong coupling assumptions on cri- 
ticality. The model was simulated in [6]; with 
the available statistics the results did not show a 
clear or conclusive signal for a critical behaviour 
for Kc(0) < |. Rather they indicate that there 
might be no critical point but just a cross-over 
like behaviour. So either there is indeed no criti- 
cal point or the available statistics is by far not 
sufficient. 

2. PARTITION FUNCTION ZEROS 

An alternative approach to investigate the 
phase structure is to determine the partition func- 
tion of the Schwinger model for fixed /3 and k G C 
and to study the finite size scaling behaviour of 
the closest Lee-Yang zeroes [9,10]. The partition 
function is given by 

Zs,a{k) = y"rf//((7)detMA(K,(7)e-^°('''^) 

= ZG,A(detMA(K,;7)) , (2) 

which is a polynomial of degree 2|A| in k. For 
a finite system the zeroes in Zs^a(k) are strictly 
complex and the partition function and the free 
energy are analytic in a non- vanishing neighbour- 
hood of the real axis. As criticality is approached 
(in the thermodynamic limit and for k Kc) the 
Lee- Yang zeroes pinch the real k axis, inducing 
a phase transition. The finite size scaling beha- 
viour of these zeroes, in particular those closest 
to the real axis, is related to the thermodynamic 
critical exponents [11,10]. 

There are three ways to determine (_D) = 



(det M\(k, U)) for complex k. 

• One obtains the coefficients of the poly- 
nomial expressed in terms of moments of 
M^^, by a direct numerical simulation. Ho- 
wever, since coefficients may not decay suf- 
ficiently fast with n, we may have to find all 
coefficients. This approach is closest related 
to the, elsewhere very successful, method of 
analytic extrapolation via histograms [10]. 

• One may calculate values (_D) with suffi- 
cient precision for various real k, perform 
a polynomial fit and use this fit to analyti- 
cally continue to complex k. In a test of this 
method for the analytically solvable free fer- 
mion theory it turned out to be very unsta- 
ble and could not be used to identify the 
closest zeroes. 

• Finally the hard but direct way: Calcu- 
late (_D) for K G C and scan certain regi- 
ons of the complex plane for the closest ze- 
roes. This method is applicable for small 
lattices only, since one has to determine the 
determinant for each gauge field configura- 
tion explicitly worked. On the other hand 
it does not involve any extrapolation and 
works for arbitrary number of fiavours. 

We used the last mentioned method in our analy- 
sis; since it is very CPU-time consuming we were 
restricted to relatively small lattices of the size 
|A| = 2^, 4^, 8^. Due to the smallness of the lat- 
tices a conclusive statement on the existence of a 
critical point may not possible. However, compa- 
ring the finites size scaling behaviour of the zeroes 
with that of the analytically known free fermion 
theory does give a hint what to expect. 

In fig. 1 we show the positions of the closest 
zeroes for the three values of /3 = 0, 1 and 5. 
A direct comparison with the values at /3 = cxd 
is not possible, since there the absolute positions 
depend on the boundary conditions. (For finite 
/3 and U(l) gauge group the b. c. for the fermi- 
ons are irrelevant and averaged out anyhow). In 
fig. 2 we compare the scaling behaviour of the 
imaginary part of the closest zeroes with that for 
the free fermion theory. We find no significant 
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Figure 1. Position of the closest zero in the com- 
plex K plane for lattice sizes L = 2, 4, and 8 and 
gauge coupling /3 = 0, 1, and 5. 



difference and observe a scaling compatible with 
Im(Ko) cx l/L. 

All data are consistent with the assumption 
that there is a critical Kg for these gauge coup- 
lings. This is particularly interesting for /3 = 
where the eight-vertex data [6] suggested that 
there is no critical point at finite k(0). On the 
other hand, much larger lattices were studied 
there. Further studies with enhanced statistics 
may clarify this discrepancy. 

The phase structure extrapolated from the sca- 
ling behaviour of the zeroes is thus in agreement 
with the assumption that there is a critical line 
starting at Kc{oo) = ^ continuing through all va- 
lues of f] and ending at a Kc(0) < |. 
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Figure 2. The imaginary part of the closest Lee- 
Yang zero Im(Ko) as function of the lattice size 
L. (a) A In-ln plot, comparing values at /3 = 0, 1, 
and 5 with those at cxd. (b) Scaling plot vs. l/L, 
the scaling behaviour expected for /3 = cxd. 
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